!
!  phaphot     Aperture photometry
!
!  Copyright © 2010 - 2011 F.Hroch (hroch@physics.muni.cz)
!
!  This file is part of Munipack.
!
!  Munipack is free software: you can redistribute it and/or modify
!  it under the terms of the GNU General Public License as published by
!  the Free Software Foundation, either version 3 of the License, or
!  (at your option) any later version.
!  
!  Munipack is distributed in the hope that it will be useful,
!  but WITHOUT ANY WARRANTY; without even the implied warranty of
!  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!  GNU General Public License for more details.
!  
!  You should have received a copy of the GNU General Public License
!  along with Munipack.  If not, see <http://www.gnu.org/licenses/>.


program phaphot

  use ifits

  implicit none

!  real, dimension(12) :: apertures
!  real, dimension(2) :: ring
  character(len=1024) :: sb
!  integer :: i

!  forall( i = 1:size(apertures) )
!     apertures(i) = 2.0*exp(0.1567271184*((i-1)*1.570796327))
!     ! logarithmic spiral, pitch 8.9073 deg
!  end forall
!  ring = (/20.0, 30.0/)

  do
     read(*,*,end=20) sb

     call faphot(sb)!,apertures,ring)

  end do

20 continue

!contains


end program phaphot

